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ABSTRACT 

It was recently found [1,2] that in sulphur-induced nuclear collisions at 
200 A GeV the observed strange hadron abundances can be explained within 
a thermodynamic model where baryons and mesons separately are in a state 
of relative chemical equilibrium, with overall strangeness being slightly un- 
dersaturated, but distributed among the strange hadron channels according 
to relative chemical equilibrium with a vanishing strange quark chemical po- 
tential. These studies were made either in a narrow central rapidity window 
[1] or (in the case of S-S collisions [2]) also using the rapidity-integrated total 
multiplicities. 

We develop a consistent thermodynamic formulation of the concept of 
relative chemical equilibrium and show how to introduce into the partition 
function deviations from absolute chemical equilibrium, e. g. an undersatu- 
ration of overall strangeness or the breaking of chemical equilibrium between 
mesons and baryons. We then proceed to test on the available data the hy- 
pothesis that the strange quark chemical potential vanishes everywhere, and 
that the rapidity distributions of all the observed hadrons can be explained 
in terms of one common, rapidity-dependent function /U q (ry) for the baryon 
chemical potential only. The aim of this study is to shed light on the ob- 
served strong rapidity dependence of the strange baryon ratios in the NA36 
experiment [3]. 



*Work supported by GSI, DFG and BMFT. 
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INTRODUCTION 



In the meantime it is almost commonly agreed that fireballs resulting from collisions 
between heavy ions probably live long enough to establish a state of local thermal and at 
least relative [1,2] chemical equilibrium. But according to our knowledge no thermodynamic 
description has been developed yet which gives a quantitative explanation for the different 
shapes of the rapidity distributions. 

We introduce a simple hydrodynamical model of fireballs resulting from S+S collisions at 
200 A GeV and allow as a new feature a dependence of the chemical potential on the space- 
time rapidity. We show that this simple extension increases dramatically the agreement 
with experimental spectra. 

We use the complete set of rapidity distributions for n~, K^, K° s , A, A and p-p from 
200AGeV/c S+S collisions obtained by the NA35 collaboration [4]. Where applicable we 
back-extrapolate the published data to the experimental acceptance windows (by inverting 
the algorithm given by NA35 [4]) since our model is able to take them into account explicitly. 
In the figures below the index "e" is short for "within experimental acceptance" . 

RELATIVE CHEMICAL EQUILIBRIUM 

J. Rafelski et al. [1] suggested a simple method to test experimental data for possible 
signals of a QGP in the early stages of heavy ion collisions based on a chemical analysis 
of the final hadronic composition. The underlying concept rests on the assumption of local 
thermal and absolute chemical equilibrium with respect to up- and down-quarks, at least 
in a region near central rapidity. Absolute chemical equilibrium of strange particles is 
not expected due to the rather high mass of ss-pairs. Instead one assumes the weaker 
condition of relative chemical equilibrium: The strange phase space is not fully saturated, 
but strangeness has been distributed among the available strange hadron channels according 
to equilibrium fugacities, maximizing the entropy. A convenient way of parametrizing this 
is by introducing a common saturation factor < 7 S < 1 for both strange quarks and anti- 
quarks. The abundance of hadrons is then regulated by fugacities X k given by a product 
of valence quark fugacities multiplied by the factor 7<!, Sfc ', where \sk\ counts the number of 
strange quarks plus anti-quarks in the respective hadron species. 

In [1] the factor 7 S was introduced heuristically; here we present a rigorous thermody- 
namic formulation of the concept of relative chemical equilibrium. 

In a statistical description of an ideal gas mixture containing M different components, 
the one particle distribution functions f k = f k (x,p), k = 1 . . . M contain all of the system's 
information. Given sufficient time, this distribution function will evolve to the local equilib- 
rium solution, whose form can be very efficiently derived by maximizing the entropy subject 
to constraints from conservation laws. For thermal and absolute chemical equilibrium the 
only constraints are due to the absolutely conserved quantities energy, baryon number and 
(on the time scale of nuclear reactions) strangeness. (Charge conservation in practice doesn't 
play a role in nuclear collisions.) However, if different chemical processes occur at different, 
well separated time scales, a state of relative chemical equilibrium may arise at intermediate 
time scales where the fast processes have already equilibrated but the slow ones are still 
far off equilibrium. This intermediate state can be characterized by additional constraints 
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on the chemical composition which, when implemented into the procedure of entropy max- 
imization, yield the local equilibrium distribution for relative chemical equilibrium. 

Let us consider a small fluid cell at freeze-out in its own rest system. Its entropy content 
can be written as 

s = -J2 I duj{f k in f k + e k (i- e k f k ) in (i - e k f k )} (i) 
k J 

where the integral is over conventional phase space / du = J AV d?x J -^s, and 6 k = +1(— 1) 
applies to baryons (mesons). The conservation of energy, light (q) and strange (s) quark 
flavours in strong interactions leads to the constraints 

(E) = I duJ2E k f k , (Q) = f duJ2qkfk, ~(S) = f duY^Skfk. (2) 
J k J k J k 

Here (E) denotes the cell's average thermal energy, (Q) stands for its average net light quark 
number (we neglect the small u-d mass difference), and (S) for its average amount of net 
strangeness. On the right sides q k and s k count the number of light resp. strange valence 
quarks in a hadron of species k; note that anti-quarks contribute with negative signs. 

The concept of a partially saturated strange phase space can be implemented by the 
additional condition 

(\S\)= f duJ2\sk\fk. (3) 
J k 

Since the amounts (S) and (\S\) of overall strangeness can be specified independently, this 
allows for a description of both over- and under-saturated strange phase space. We can 
further parametrize 

(\s\} = T,^(\s\)k\ (4) 

k 

where (l^l)^ 01 denotes the contribution of the hadron species k to overall strangeness in 
absolute chemical equilibrium. The factors represent the degree of population of the 
strange phase space in each hadronic sector, which may differ from their equilibrium values 
= 1; the range < < 1 corresponds to strangeness suppression. The parametrization 
(||) explicitly takes into account that quarks are bound into hadrons at freeze-out and allows 
for a suppression in dependence on their valence quark content. 

The most probable distribution function f k can be found by extremizing the functional 
<S[/fc], subject to the constraints (0, This can be done by means of Lagrangean multipliers 
P, u q , u s , a and leads to 

^ k ~ e PE k +v q q k +v s s k +a\s k \ ZjZ Q k ' ^ 

Inserting (|5|) into the entropy formula (|l]) we obtain 

S = P(E) + ^(v q q k + v a s k )(N k ) + l*E + Y l <r\sk\(Nk) , (6) 
k k 
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where we define the average particle number of species k in the considered volume element 
(N k ) = J dujf k and 

In 5 = ^ 6 k f do; In (l + e k e- WEk+i ^ qk+ " sSk)+CTlSkl) ) . (7) 
k J 

The limit a — > corresponds to neglecting the additional condition (|3|) in the extremization 
procedure. In this case we recover from fl7|) the expression for the grand canonical partition 
function in absolute chemical equilibrium. 

It can be easily shown that InS generates all thermodynamic quantities (ac k = (v q qk + 
u s s k ) + a\s k \): 

(E), ~lnB a =W, Tj-lnH =p. (8) 



- — In 3 
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Hence it is a thermodynamic potential. As usual we identify /3 = 4 and a*. = — ^r. The 
chemical potentials in absolute chemical equilibrium 

^Eq _ _T^ qk + z/sSfc ) (g) 

are connected with /i^ via 

/i k = fi k Q - Ta\s k \ = fi k Q + a\s k \ , (10) 

where we defined a = —Ta. 

In Boltzmann approximation we can determine the suppression factors by comparing 
(0) and (ED) as = e^ CT ' Sfe '- ) . They are directly related to the saturation factor 7 S via ^ k = "fl Sk ^. 

0< 7 M = efe = (e^) H <l. (11) 

An overall suppression of the strange quark phase space has no influence on the thermody- 
namic relations (§); it only enters into the relation ([!]) for the entropy. If decomposed into 
particle specific contributions S = J2k^k we immediately see that each component receives 
an additional contribution —\s k \ln / y s (N k ). 

The starting point of our formalism were the suppression factors weighting each hadron 
species separately. As a consequence of the entropy maximization, however, we automatically 
arrived at the concept of relative chemical equilibrium with one common factor j s . This 
justifies a posteriori the parametrization of Rafelski et al. [1]. 

The procedure is easily generalized to allow for a meson-baryon non-equilibrium as sug- 
gested in the last of Refs. [1]. This might be necessary if processes of the type M+M «-> 
B+B are slow. Let us assume that baryons are in absolute chemical equilibrium and mesons 
are not. This requires the inclusion of the additional constraint 

(M) = J2 \{i-e h ) I dwj k , (12) 

k J 

where again we can parametrize in the same spirit as (f|) 

W = El(i-^)4 m W Eq . (is) 

k 
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Combining (|3]) and (|TJ) we find the relative equilibrium distribution 

fu = il Skl 4 {1 ~ 9k) A Eq , (14) 

where in Boltzmann approximation £™ and 7 m are related by £™ = 7m 1 ° k ■ The case where 
also the baryons are suppressed or enhanced relative to their equilibrium abundance can be 
treated in an analogous way. 



CALCULATING THE SPECTRA 

For the momentum spectra of hadrons emitted from the fireball at freeze-out we use the 
simple hydrodynamical model developed by E. Schnedermann et al. [5]. It is based on the 
idea of a thermal distribution with one freeze-out temperature Tf characterizing the instant 
at which fluid cells loose their thermal contact to the fireball. Experiments find similar 
exponential slopes in all m_i_-spectra, independent of the rapidity window [4]. This supports 
the interpretation of freeze-out at one common and constant temperature. 

The general expression for the invariant momentum distribution is [6] 

E d ^r = 7 ^T-J f(x,p)p"d 3 a u , (15) 



dp 3 (2vr) 3 Ja { 

where g is the spin degeneracy of the hadrons under consideration and <7f is the 3-dimensional 
freeze-out hypersurface in space-time, along which the condition T = Tf is fulfilled. We 
choose a cylindrical coordinate system (r, <j), z) which best reflects the global symmetries of 
the expanding fireball and obtain for the normal vector of the freeze-out surface 

d 3 a u = ( 1, — ^ cos0, — ^ sin0, — ) rdrd(pdz. (16) 
\ Or Or Oz J 

In an azimuthally symmetric geometry of this kind it is practical to decompose the velocity 
field in the following way [7]: 

u u = (cosh p cosh rj, e± sinh p, cosh p sinh rj) . (17) 

Here rj = t](t,r,z) is the longitudinal flow rapidity by which each volume element on the 
z-axis moves relative to the center of mass, and p = p(t, r, z) is the rapidity corresponding 
to the transverse flow of the volume element at position (r, z) as seen from a reference point 
at z on the beam axis moving with the local flow velocity there. The 4-momentum can be 



parametrized in terms of rapidity y and transverse mass m± = \jp\ + Wq as 

p v = (m± cosh y, p± cos P , p± sin (fi p , m± sinh y) . (18) 
After suitably orienting the coordinate system we obtain in Boltzmann-approximation: 
d 2 N 



dydm\ 



-^-m ± ( dzrdre^ e- &cosh{ri - y) 

47T Jt { 



X 



coshy - sinhy^ j I (a) - > ( 19 ) 
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where a = (m±/T) cosh p and a = (p±/T) sinhp. The 2-dimensional integration extends 
along the curve determined by all points (r, z) which at time tf(r, z) satisfy T(t{ (r, z),r, z) = 

For further evaluation we introduce new variables and transform (t, r, z) into (r, r, Y) 
with the longitudinal proper time r = \Jt 2 — z 2 and the space-time rapidity Y = 
0.5 In \(t + z)/(t — z)]. Numerical simulations of the space-time evolution of the hot zone 
in ultrarelativistic nuclear collisions show an almost linear velocity profile in beam direc- 
tion [8]. Hence a boost invariant scenario, if restricted to the fireball's finite extension, 
is a good approximation for the longitudinal fluid dynamics and leads to Y = 77, i.e. the 
identity of space-time rapidity Y and fluid rapidity 77. The partial derivatives |p and ^ L 
are easily evaluated if we parametrize the freeze-out hypersurface by a proper time surface 



Tf(r): tf(r,z) = yT 2 (r) + z 2 . In doing so we simultaneously express the invariance of the 
decoupling process against longitudinal boosts and obtain for the rapidity distribution using 
z = Tf(r) sinhr?: 



dN g f m t , 2 pi J fRf 
dy Atc 



/ dm\ \ drj drT { (r)re^ T e- &co8h ( y -ri (20) 

Jm}° J-rj Jo 

( dr i(r) \ 

x I m ± cosh(y -r/)I (a) - p±_ — Ii(«) 1 . 

In this expression m 1 ^ , nr^ denote the experimental limits in which the spectrum was mea- 
sured, the freeze-out radius is called R{ and the longitudinal extent of the fireball is fixed via 
the finite interval — fj . . . 77. We will later use the approximation where we assume transver- 
sally instantaneous freeze-out ( 8T gy = 0) and replace the transverse flow velocity profile by 
its radial average. In this case eq. (|20|) reduces to the simpler form 



dN _ g TfRf 
dy 8 7r 



2 „ m h 



\ dml m± I (o;) / d V cosh(?/ - 77) e^ /T e - &cosh ^ . (21) 

Jm'f J— fj 



THE PARAMETRIZATION OF THE CHEMICAL POTENTIAL 

Spectra of hadrons k containing up or down valence quarks and Sk strange quarks are 
fully characterized, besides the two flow components, by the temperature T = 1/(3 and the 
fugacity = 7j, Sfc ' e^ MsSfe e^^ qk (see the section on relative chemical equilibrium). 

J. Sollfrank et al. [2] found in a chemical analysis of S+S collisions, using rapidity inte- 
grated total multiplicities, a vanishing chemical potential fi s ~ of the strange quarks and 
an almost complete saturation 7 S w 1 of the phase space of strange particles. It is doubtful 
that these parameters can be established in a purely hadronic environment, therefore we 
adopt as a working hypothesis the formation of a QGP in the early stages of the collision. 
This leads to /x s = independent of T, fi q and also of r and z (resp. 77). We further assume 
that these values are maintained during the process of hadronisation. (For a detailed dis- 
cussion of this assumption see the last reference [1].) Our objective is to completely fix the 
baryon chemical potential by using only the hyperon spectra, and then test whether other 
measured hadron distributions can be reproduced in shape as well as normalization. 
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FIG. 1. Rapidity distribution of A- and A-hyperons. The theoretical curve with our quartic 
parametrized chemical potential is fitted to data of NA35 S+S 200 AGeV [4]. 

The width of the rapidity distributions is mainly determined by the longitudinal flow 
component. It can be extracted with very little sensitivity to T and transverse flow from 
the shape of the pion rapidity spectrum [5]. In this way we obtain a maximum fluid rapidity 
(see (|21~D) of fj = 1.75. In a first attempt we neglect the transverse flow (p = 0) and choose 
T = 200 ± 20 MeV, consistent with the slope of transverse mass distributions. We are aware 
of the fact that this value is too high to be interpreted as the real freeze-out temperature 
directly, but in a hadron gas a temperature of this order seems to be necessary to ensure 
strangeness neutrality [1,2,9]. 

Both protons and A-hyperons possess a quite broad rapidity distribution and peak in 
the fragmentation regions. This suggests a concave parametrization for n q (r]). We checked 
several possibilities and found that a quartic ansatz of the form 

^(rj) = + Ar] 4 )q k , fj° = 59 ± 11 MeV , A = 11.5 ± 0.5 MeV (22) 

together with fj = 1.75 is able to describe both A- and A-distributions simultaneously 
(Fig. [I] a-b). The central data points of the A's are barely reached only, but they have 
the largest error bars. We used the fact that both spectra are measured in the same p±- 
window, whereby the integrals over the boostangle r] cancel in the ratio of the distributions 
at midrapidity. From the central normalization of the A-distribution we also determine 
TfRf = 20±4fm 3 . The contribution from S° decays to the A spectrum is taken into account 
in the theory by an additional degeneracy factor of 2. 



THE RESULTS 



Integrating fi(rj) over the fluid (or space-time) rapidity t], we determine an average baryo- 
chemical potential of /2b = 3/2 q = (240 ± 40) MeV. This is in good agreement with the 
findings of J. Sollfrank et al. [2], /ib = (264 ± 72) MeV who used a chemical analysis 
of the 47r-integrated multiplicities. But the most important result is that we are able to 
also reproduce all other available hadron spectra (Fig. |2| a-d). Not only the shape of the 
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distributions coincides with the experiment, but also (except for the pions which will be 
discussed separately below) the normalization is reproduced correctly without adjustments. 
(Due to error propagation some of our kaon curves may, however, have up to 35% uncertainty 
in normalization.) This is a dramatic improvement over the results obtained in Ref. [5] with 
rapidity independent chemical potentials. 

The calculated proton distribution (Fig. |2| d) is clearly narrower than the experimental 
result, moreover the normalization is a little too low. The first feature is expected because 
protons are unlikely to be completely thermalized, especially in the extreme kinematic do- 
mains, due to nuclear transparency in particular in the transverse nuclear surface. On the 
other hand a slight undersaturation of the overall strangeness by a factor 7 S = 0.75 fully 
corrects the normalization problem (Fig. |2] e). Since all the other hadrons carry only one s- 
or s-quark, this leaves the spectra unchanged, if the factor Tf R| is increased correspondingly 
by a factor 4/3. 

We conclude from this agreement with experiment that a thermodynamic description 
will work well in a rather wide rapidity interval if dependencies of the thermodynamic 
parameters on rapidity are taken into consideration. The interesting observation here is 
that it is apparently sufficient to consider a rapidity dependence of the baryon chemical 
potential, keeping T, fi s = and 7 S fixed. On the other hand this means that one has 
to be careful when applying thermodynamic models to data which are not restricted to a 
narrow and central kinematic domain, because then this rapidity dependence is integrated 
over with different weights for different particle species. As a further result the strong 
dependence of strange particle ratios on the rapidity window which was observed by NA36 
in S+Pb collisions [3] can be qualitatively explained by our model in a very natural way. 

To test the important condition of strangeness neutrality we define 



e(ry) 



(23) 



as a measure for the local deviation. Fig. [3] shows the result plotted in the relevant range of 77 
for two different cuts in the hadron mass spectrum. We recognize that strangeness neutrality 
can not be guaranteed locally, though e(rj) is flat and almost zero in —1.2 < rj < 1.2. On 
the other hand the clear experimental differences between the rapidity densities of different 
hadrons can be seen as an indication that a violation might indeed be realized locally. In a 
global sense, i.e. integrated over rapidity 7], we obtain 



7s 


1.0 


0.75 


m cut [GeV] 


1.8 


2.0 


1.8 


2.0 


^global 


0.03 


-0.16 


0.09 


-0.09 



If we take into consideration the uncertainties in temperature, chemical potential and a 
possible asymmetry in the data tables between meson and baryon resonances of high rest 
masses (due to different efficiencies for their identification), we conclude that we have no 
problem in ensuring at least global strangeness neutrality. 
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FIG. 2. Rapidity distribution 
of kaons and protons in compar- 
ison with data of NA35 for S+S 
at 200AMeV [4]. The absolute 
normalizations are not fitted, but 
calculated using T = 200 MeV, 
fj = 1.75 and the quartic 
parametrization for the chemical 
potential. Figure e) represents 
the proton distribution assum- 
ing an undersaturation of overall 
strangeness 7 S = 0.75. 
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FIG. 3. Deviation from local 
strangeness neutrality. The quan- 
tity e (23) is plotted func- 
tion of space-time rapidity r\ for 
T = 200 MeV, 7 S = 1 and the 
parametrization (22) for the chem- 
ical potential. The dashed curve 
is calculated using a mass cut in 
the hadronic resonance spectrum 
at m cut = 1.8 GeV, the solid one 
represents m cu t = 2.0 GeV. 



-1 1 T| 

TRANSVERSE FLOW 

Despite this very good and consistent agreement of our model with experiment this 
comparison has one substantial shortcoming: the assumed temperature is really too high 
for consistency of the model. The spatial overlap of hadrons at T = 200 MeV is quite large 
already, and indeed lattice QCD simulations now firmly predict a QGP phase transition 
already at T = 150 MeV [10]. The assumption of a system of non-interacting particles 
freezing out at T = 200 MeV is thus more than dubious. These were the decisive reasons for 
the authors of [2,5] to conclude that the flat slope of the observed mj_-spectra is partially 
caused by a blueshift due to collective transverse flow. Already p ~ 0.3 reduces the true 
temperature from T ~ 200 MeV to T ~ 150 MeV, and we will use this average value in the 
following considerations. 

We again determine and A as described above from the hyperon spectra, now leading 
to A = 8.0 ± 0.5 MeV and ~ 43 MeV. The lower temperature reduces the thermal 
smearing whereby all the theoretical rapidity spectra become somewhat narrower. The 
calculated proton distribution now fits the experimental data (except for the disagreement 
near the projectile and target rapidities discussed already above), even with 7 S = 1 (Fig. [| c). 

With kaons, however, we now have a serious problem: Although the shape is correct 
for all three charge states, the theoretical normalization exceeds the data by a common 
factor of approximately 2.7 (Fig. |] d-e). As a consequence the theoretical spectra also 
strongly violate overall strangeness neutrality. However, there is a very simple remedy for 
this. Following the ideas first laid out in Ref. [Id] and worked out in Ref. [lc], one can 
assume a breaking of the relative chemical equilibrium between the mesons and baryons 
by introducing a meson/baryon suppression factor 7 m < 1. It is quite a surprise that by 
choosing 7 m = 0.37 we can simultaneously correct the normalization of all kaon spectra 
(Fig. H) and recover global strangeness neutrality: 
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tions of kaons and protons in com- 
parison with data of NA35 S+S 
at 200AMeV [4]. The absolute 
normalizations in the cases c)-e) 
are not fitted, but calculated us- 
ing T = 150 MeV, fj = 1.75 and 
the quartic parametrization for the 
chemical potential. Please note 
that 7s = 7m = 1 was assumed. 
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FIG. 5. Rapidity distribution of kaons in comparison with data of NA35 S+S 200AGeV [4]. 
The theoretical curve is calculated using T = 150 MeV, fj = 1.75 and 7 m = 0.37. 

EXCESS ENTROPY 

In both analyses (with and without transverse flow) we used the pion rapidity spectra 
to extract the maximum flow rapidity, but did not check their normalization. It turns out 
that in both cases we seriously underpredict the pion multiplicity (typically by a factor 2-3). 
While in the case with transverse flow and 7 S = 7 m = 1 we are not too much off (see fig. ^ a), 
after introducing the meson suppression factor 7 m = 0.37 (thereby recovering the correct 
kaon normalization and strangeness neutrality) we again obtain by a factor 2-3 too few pi- 
ons. This feature in the data (too many pions compared to a hadron gas in relative chemical 
equilibrium) was observed before [1,2,11] and is confirmed by the completely different ex- 
perimental analysis reported by M. Gazdzicki at "Strangeness '95" [12]. It is interpreted as 
an indication for excess entropy in the data which might originate from the gluonic entropy 
in the QGP [1]. Why this extra entropy shows up mainly in the form of pions is not yet 
clear. Their nature as nearly massless chiral goldstone bosons may play a role here, but no 
dynamical explanation has been given so far. 

SUMMARY 

We showed that a rather simple thermodynamic model with a suitably parametrized 
rapidity dependent baryochemical potential can reproduce the measured distributions of 
strange particles resulting from S+S collisions at 200 A GeV. Consistency was achieved with 
the shape as well as with the normalization of the rapidity spectra. Our phenomenologi- 
cally motivated ansatz was not able to guarantee strangeness neutrality locally, but globally 
strangeness was conserved. One further result of our model was a qualitative explanation of 
the strong dependence of strange particle ratios on the chosen kinematic window. 

Resonance decays have not been taken into consideration so far, but will be considered 
in the future. They are expected to mostly affect the pion spectra, but also kaons could 
be influenced in a non- negligible way. We have conscisously kept /i s = and 7 S , 7 m as 
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FIG. 6. Rapidity distribution of pions. The theoretical curve with our quartic parametrized 
chemical potential and a transverse flow of /?t = tanh p = 0.3 is compared to data of NA35 S+S 
200AGeV [4]. Figure b) is based on a meson baryon suppression factor j m = 0.37. 

rapidity- independent, because otherwise a unique determination of such a large set of free 
parameter functions would be impossible to fix from the available data. 
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